Giant halo in relativistic and non-relativistic approaches 
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Abstract 


The phenomena of giant halo and halo of neutron-rich even-Ca isotopes are investigated and 
compared in the framework of the relativistic continuum Hartree-Bogoliubov (RCHB) and non- 
relativistic Skyrme Hartree-Fock-Bogoliubov (HFB) calculations. With two parameter sets for 
each of the RCHB and the Skyrme HFB calculations, it is found that although halo phenomena 
exist for Ca isotopes near neutron drip line in both calculations, the halo of the Skyrme HFB 
calculations starts at a more neutron-rich nucleus than that of the RCHB calculations, and the 
RCHB calculations have larger neutron root-mean-square (rms) radii systematically in N > 40 
than those of the Skyrme HFB calculations. The former difference comes from difference in shell 
structure. The reasons for the latter can be partly explained by the neutron 35} /2 orbit, which 
causes more than 50 % of the difference among the four calculations for neutron rms radii at ©°Ca. 


PACS numbers: 21.10.Dr, 21.10.Gv, 21.60.Jz, 27.40.42, 27.50.+e 


I. INTRODUCTION 


The study of exotic nuclei with extreme N/Z ratios has attracted worldwide attention 
since the first discovery of halo in "Li [1]. Usually for exotic nuclei, their Fermi surfaces are 
very close to the continuum threshold, and the valence nucleons could be easily scattered 
to the continuum states due to the pairing correlations [2, 3]. Thus, theories which can 
properly handle the pairing and continuum states are needed to describe the properties of 
exotic nuclei and to understand the interference of the continuum-energy states to bound 
many-body systems and effects of large neutron excess. 

It is well known that the BCS approximation for the pairing correlations suitable for 
stable nuclei cannot give correct wave functions near the drip lines, see e.g. [4]. Taking 
into account the pairing correlations by the Bogoliubov transformation, relativistic and/or 
non-relativistic mean field approaches have been extensively used to describe and predict 
halo phenomena in exotic nuclei. 

The status of the studies of the halo using the relativistic continuum Hartree-Bogoliubov 
(RCHB) or the non-relativistic Hartree-Fock-Bogoliubov (HFB) method include: the pre- 
dictions of halo in the Ne [5, 6], Na [6, 7], Ca [8] and Zr [9] isotopes near the neutron drip 
line. More details can be found in recent review papers, e.g., [10]. As there are more than 
two particles in the weakly bound or the positive-energy region in Ca and Zr isotopes near 
the neutron drip line, the halo phenomena for Ca and Zr are addressed as giant halo [8, 9]. 
The proton-magic nuclei from O and Pb nuclei are also systematically surveyed, and it is 
found that the halos signature in Zr is the clearest in Ref. [11]. There is no evidence for halo 
in Ni, Sn and Pb isotopic chains, similarly as in Ref. [12]. 

In the non-relativistic HFB approaches, Skyrme calculations with the parameter set SLy4 
predicts the halo in Sn and Ni [13]. It is also claimed that the pairing gaps have an effect 
to reduce the halo at the neutron drip line [13, 14, 15}. 

For comparison with the halos predicted in relativistic approaches |8, 9, 11], it is necessary 
to investigate the Ca and Zr in non-relativistic approaches, which is very important but still 
missing in the literature. Therefore the present paper is dedicated to investigate the neutron- 
rich even Ca isotopes up to the neutron drip line. The results from both relativistic and 
non-relativistic mean field approaches taking into account the pairing correlations by the 
Bogoliubov transformation will be presented and compared in detail. 

The numerical details used in both calculations are explained in Sec. II, and then the 
neutron-number dependence of the energy and the rms radii is examined in Sec. III. The 
properties of single particles are discussed in Sec. IV in order to understand the halo struc- 
ture. A brief summary is given in Section V. 


Il. METHOD OF CALCULATION 


The detailed formulations of the RCHB method and the HFB method can be found 
respectively in Refs. [2, 3] and [16, 17]. Both calculations are based on the coordinate 
representation and use a box size 20 fm with the spherical symmetry assumed and a mesh 
size 0.1 fm. The quasiparticle states are obtained up to 120 MeV (165-190 radial wave 
functions in ®°Ca), and all of these states are used for calculating potentials. The maximum 
angular momentum of the quasiparticles jmax is = (RCHB) and 2 (Skyrme HFB). The 
parameter sets NL-SH [18] and PK1 [19] ( SkM* [20] and SLy4 [21] ) are used for the RCHB 
( Skyrme HFB ) calculations. For the pairing part, a surface-type delta interaction [4, 9] 
is used with the density parameter pọ = 0.152 fm~? and 0.160 fm~? for the RCHB and the 
Skyrme HFB calculations, respectively. 
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FIG. 1: Binding energy per nucleon B/A of even Ca in RCHB calculations with NLSH and PK1 
and the Skyrme HFB calculations with SkM* and SLy4 compared with the experimental data [23] 
at N = 30 and 32. For the details, see the text. 
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FIG. 2: Two-neutron separation energy Sən of even Ca in RCHB calculations with NLSH and 
PK1 and the Skyrme HFB calculations with SkM* and SLy4 compared with the experimental data 
[23]. Note that Sən of SLy4 are very close to those of PK1 in 42 < N < 48. 


The strength of the pairing interaction in mean-field calculations is normally deter- 
mined so as to reproduce the pairing gaps obtained from odd-even mass differences (see 
e.g. Ref. [22]). While for nuclei close to the neutron drip line, the pairing energy obtained 
using the Gogny interaction in the pairing channel is adopted to fix the strength of surface- 
type delta pairing interaction due to the missing of experimental data. The Vo is —325 
MeV fm? in the RCHB calculations, and that of the Skyrme HFB calculations is —365 
(—330) MeV fm? in the SLy4 (SkM*) calculations. The average pairing evergy of the PK1 
and NL-SH calculaitons —8.2 MeV at °°Ca is used for determining the Vo of the Skyrme-HFB 
calculaitons. We checked that when jmax = 8 was used in the Skyrme HFB calculations, 
the rms radii did not change with the re-adjusted Vo from the pairing energy. 


II. NEUTRON-NUMBER DEPENDENCE OF THE HALO STRUCTURE 


Figures 1 and 2 show the binding energy per nucleon B/A and the two-neutron separation 
energy Sən in RCHB calculations with NLSH and PK1 and the Skyrme HFB calculations 
with SkM* and SLy4, respectively. The most important differences between these calcu- 
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FIG. 3: The neutron and proton rms radii 7, and 7, of even-Ca isotopes in RCHB calculations 
with NLSH and PK1 and the Skyrme HFB calculations with SkM* and SLy4 in the neutron-rich 
region. 


lations are location of the two-neutron drip line and a shell gap at N = 40. The RCHB 
calculations show the two-neutron drip line at N = 50 (PK1) or 52 (NL-SH), and the pre- 
diction of the Skyrme HFB calculations is N = 46 (SLy4) or 58 (SkM*). The RCHB and the 
SLy4 calculations have clear the shell gap at N = 40 (Fig. 2), while the SkM* calculation 
does not show at all. The SkM* calculation has a small extra downward behavior of So, at 
N = 50, however, it is too small to call a shell gap. Therefore, no new shell gap is predicted 
in N > 40, in any of the calculations. The two RCHB calculations give the results rather 
close to each other both in B/A and Sn, and the two Skyrme HFB calculations have a 
difference of 0.1-0.2 MeV in B/A. The Sən of SLy4 is close to those of the SkM* calculation 
in 34 < N < 40 and then coincides with that of PK1 in 42 < N < 48. The experimental 
data are available up to N = 32 currently [23]. The RCHB calculations reproduce well 
the data of B/A at N = 30 (8.550 MeV) and 32 (8.396+0.013 MeV) and Sən at N = 32 
(9.081+0.699 MeV). For Sən at N = 30, the SLy4 calculation gives a value very close to the 
measured one (11.499+0.008 MeV). 

Figure 3 shows systematics of the rms radii of neutrons 7,, and protons rp of even-Ca from 
N = 30 to the two-neutron drip lines. The increase in the curvature of 7,, indicates the halo. 
(For attempts to define the criterion of halo or the halo size, see Refs. [7, 13, 24] and references 
therein.) For Ca, the number of nucleons in the positive-energy region N, of °~"Ca is 0.6- 
2.2, of which average is 1.7, in the NL-SH calculation [8], and the corresponding value of the 
SkM* calculation is always smaller than 0.5 (the average 0.27). As there are more than two 
nucleons in the weakly bound orbits and the positive-energy region for these nuclei and also 
halo phenomena in neighboring nuclei with incremental neutrons, the halo in Ca isotopes 
are refereed as giant halo, as in Refs. [8, 9]. The SkM* calculation has the halo from N = 52, 
and the RCHB calculations show the gradual occurrence of the giant halo. Apparently, the 
starting nucleus of the halo of SkM* corresponds to the extra lowering of the Sən. The SLy4 
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FIG. 4: The difference between neutron and proton rms radii 7, — Tp in RCHB calculations with 
NLSH and PK1 and the Skyrme HFB calculations with SkM* and SLy4 for even-Ca isotopes with 
N > 40. 
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FIG. 5: The difference of neutron rms radii between neighboring even-Ca isotopes with N > 40 
Tn( N) —7,(N — 2) in RCHB calculations with NLSH and PK1 and the Skyrme HFB calculations 
with SkM* and SLy4 in the neutron-rich region. The dashed line was obtained from 7,(N) = 
1.139N 1⁄3 fm [8]. 


calculation does not have the halo, because the particle-stable region ends at N = 46. The 
large difference between 7,, and 7, shown in Fig. 3 is identified with neutron skin in the 
region with no halo. 7, — 7, is displayed in Fig. 4, in which the differences between the 
parameter sets are clear; the curve of SLy4 is parallel to that of SkM*, and PK1 has a larger 
curvature than the other parameter sets in N < 50. An important aspect of the halo is 
how 7,, changes when two more neutrons are added, thus it is worth showing the difference 
Trn(N) —7,(N — 2) (Fig. 5). When there is no halo, two additional neutrons increase 7,, by 
0.05 fm, and the increase rises to 0.11 fm in the halo nuclei. In all of Figs. 3-5, the difference 
in 7, is clear between the RCHB and the Skyrme HFB calculations. It is to be noted that 
Tp is not constant, reflecting the self-consistency of the neutrons and protons interaction. 


IV. SINGLE-PARTICLE STRUCTURE 


To understand the halo phenomena and the different predictions in RCHB and Skyrme 
HFB calculations, the nucleus Ca is taken as an example. Figure 6 depicts the single 
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FIG. 6: The neutron single particle level e°? for ®°Ca in RCHB calculations with NLSH and PK1 
and the Skyrme HFB calculations with SkM* and SLy4. 


particle level e¢?" in RCHB calculations with NLSH and PK1 and the Skyrme HFB calcula- 
tions with SkM* and SLy4 for neutrons in Ca, i.e., the diagonal elements of the mean-field 
Hamiltonian in the canonical basis. 

The order of the levels is essentially the same for all the calculation. Furthermore the 
spectra for RCHB with PK1 and NL-SH are quite similar. However, the level 151/2 in SkM* 
is apparently higher than those of the other parameter sets. From left to right, the shell gap 
become less apparent. The shell gaps at N = 2, 8, 20, 28, and 40 are clearly shown in the 
RCHB spectra, while the shell gaps of N = 28 and 40 do not appear in that of the SkM* 
calculation (see also Fig. 2). 

It has been pointed out [8] that the neutron 3s1/2 level played an important role in the 
giant halo. Thus, comparison is made in Fig. 7 for the rms radii of each single-particle 
orbit (canonical basis) rą in RCHB calculations with NLSH and PK1 and the Skyrme HFB 
calculations with SkM* and SLy4. The lines of 7, versus eṣ?" are quite similar to each other 
except for the 1s,/2 and the region around «“” = 0. This indicates the similarity of the 
potential wall except for the bottom of the potential wall. A magnification around £°® = 0 
is displayed in Fig. 7b. The most significant difference is 381/2: 7351/2 of PK1 is more than 
2 fm larger than that of SkM*. These differences are clearer in Fig. 8 which illustrates 7, as 
radii of circles. 

Table I shows e€", Fa, occupation probability v2, and RUT for the neutron 3s1/2 orbit, 
as well as 7, of Ca. Both 7351/2 and va make 7, of PK1 the largest one. Since the 
contribution of the orbit to 7, is given by 20561 j273s1/2; it is seen from Tab. I that 54 % 


(0.082 fm) of the difference 7,,(PK1)—7,(SkM"*) = 0.153 fm comes from the neutron 3812. It 
is also noted that the higher energy does not necessarily mean the larger spatial distribution. 

The difference in the neutron 31/2 orbit is shown also in Fig. 9 which illustrates the 
radial wave functions squared. The tail of PK1 has appreciably longer distribution than the 
others, and as a counter part, the amplitude of PK1 is smaller in the inner region. The 
number and locations of the nodes are the same for the four lines. The neutron density 
distributions have appreciable difference (Fig. 10) in accordance with the behavior of Fig. 9. 

Figure 11 depicts central potentials. These are the summation of the scalar and the vector 
potentials in the RCHB calculation, and the summation of the tọ and t3 terms in the particle- 
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FIG. 7: a) The rms radii of each neutron single-particle orbit in canonical basis rą as functions 
of its single-particle energy ¢¢°". b) Same as a) but for the region around ¢°" = 0. 


TABLE I: Properties (e%", Fa, occupation probability v2, and 


Zura ) of the neutron 351/2 orbit and 7, of ®Ca. 


E351 /2 T3s1/2 U3s1/2 N V3s1/27381/2 Tn 

[MeV] [fm] [fm] [fm] 
PK1 0.133 8.144 0.239 0.0846 4.363 
NL-SH 0.640 7.240 0.089 0.0280 4.314 
SLy4 1.939 6.856 0.016 0.0043 4.205 


SkM* 1.445 5.765 0.010 0.0025 4.210 


hole potential in the Skyrme HFB calculation, see Eq. (A.5a) of Ref. [16]. A significant 
difference is that the SkM* calculation has smaller depth than the other calculations, as is 
anticipated from Fig. 6, and another is a slight difference in the tail region (see the inset of 
Fig. 11). It should be noted that the tails of the RCHB calculations have longer distributions 
than the Skyrme HFB calculations, and PK1 has the longest one. Figure 12 displays the 
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where p(r) is the total nuclear density, and f,(7r) is the pairing density [16] of neutron. 
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FIG. 8: The rms radii of each neutron single-particle orbit in canonical basis 7, in RCHB calcu- 
lations with NLSH and PK1 and the Skyrme HFB calculations with SkM* and SLy4 displayed as 
radii of circles. The proton one 7, (the blue center circle in the on-line version) is included for a 
comparison. The grade of color corresponds to the order of Fa. 


For the RCHB calculations, the S = 1 component of the pairing interaction is removed 
explicitly. Again, the line of PK1 has the longest tail, therefore, the large 73.1 /2 of PK1 can 
be understood in terms of the tails of the central and the pairing potentials. 

Finally, we mention the importance of the pairing correlations. It turned out that when 
the strength of the surface-type pairing interaction was determined so as to reproduce the 
neutron pairing gaps of 474446°Ca, 1.0-1.9 MeV, obtained from the odd-even mass differ- 
ences of the experimental data [23], the pairing energy of the SkM* calculation was — (35-39) 
MeV in 40 < N < 50, and the kink of 7,, at N = 50 was completely smeared out (the circles 
in Fig. 13). It is noted that the pairing energy is more than 3 times larger than the values 
of the RCHB calculations 8-9 MeV. On the other hand, the volume-type pairing interaction 
Vod(r — r’) can reproduce the pairing energy of the RCHB calculation and the experimental 
pairing gaps simultaneously, and the 7,, has the kink (the triangles in Fig. 13). Therefore, 
the kink is strongly influenced by the pairing interaction. (See also Refs. [13, 14].) 


V. CONCLUSION 


In this paper, the phenomena of the giant halo and halo of the neutron-rich even-Ca 
isotopes have been investigated and compared in the framework of the RCHB and the 
Skyrme HFB calculations. With the two parameter sets for each of the RCHB and the 
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FIG. 9: Radial wave function squared |w,|? (canonical basis) of the neutron 3s1/2 of °°Ca. They 


are normalized as f dr |y,|? = 1. 
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FIG. 10: Neutron density distribution of ĉĉCa calculated with the four parameter sets. 


Skyrme HFB calculations, it has been found that although the halo phenomena existed for 
Ca isotopes near the neutron drip line in both calculations, the halo of the Skyrme HFB 
calculations started at a more neutron-rich nucleus than that of the RCHB calculations, and 
the RCHB calculations had larger neutron rms radii systematically in N > 40 than those 
of the Skyrme HFB calculations. The former difference comes from difference in the shell 
structure. The reasons for the latter can be partly explained by the neutron 3s1/2 orbit, 
which causes 50 % of the difference in the neutron rms radii among the four calculations at 
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FIG. 11: Central potential of the neutrons of ®°Ca. Inserted is a magnification of the surface 
region. 
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FIG. 12: Pairing potential of the neutrons of Ca. 
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FIG. 13: Ff, of the SkM* calculations with three pairing interactions: the surface-type pairing 
with the strength determined from the pairing energy of the RCHB calculations (squares), the 
surface-type with the strength determined from the pairing gaps (circles), and the volume-type 
determined from the pairing gaps (triangles). The lines of the second and the third ones are very 
close to each other. 
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